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Abstract 

The densities of Yang-Lee zeros for the Ising ferromagnet on the LxL square lattice are evaluated 
from the exact grand partition functions (L = 3 ~ 16). The properties of the density of Yang-Lee 
zeros are discussed as a function of temperature T and system size L. The three different classes 
of phase transitions for the Ising ferromagnet, first-order phase transition, second-order phase 
transition, and Yang-Lee edge singularity, are clearly distinguished by estimating the magnetic 
scaling exponent from the densities of zeros for finite-size systems. The divergence of the 
density of zeros at Yang-Lee edge in high temperatures (Yang-Lee edge singularity), which has 
been detected only by the series expansion until now for the square-lattice Ising ferromagnet, is 
obtained from the finite-size data. The identification of the orders of phase transitions in small 
systems is also discussed using the density of Yang-Lee zeros. 
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I. INTRODUCTION 

Yang and Lee proposed a mechanism for the occurrence of phase transitions in the 
thermodynamic limit and yielded an insight into the problem of the Ising ferromagnet at 
arbitrary temperature in an arbitrary nonzero external magnetic field by introducing the 
concept of the zeros of the grand partition function Z(T, H) (for fluid systems, Z(T, //)) in 
the complex magnetic- field (for fluid systems, fugacity) plane (Yang-Lee zeros). They Q] 
also formulated the celebrated circle theorem which states that Yang-Lee zeros of the Ising 
ferromagnet lie on the unit circle Xq = e %e in the complex x = e~ 2H / kBT plane. Below and at 
the critical temperature T c , Yang-Lee zeros cut the positive real axis at the ferromagnetic 
transition point x c = 1 (H c = 0) in the thermodynamic limit. The spontaneous magnetiza- 
tion mo is determined by the density of Yang-Lee zeros g(9) on the positive real axis, i.e., 
m = 2-Kg(6 = 0). However, above the critical temperature, Yang-Lee zeros do not cut the 
positive real axis in the thermodynamic limit. There is a gap in the distribution of zeros 
around the positive real axis, that is, g(6) = for \9\ < 9 e . Within this gap, the free energy 
is analytic and there is no phase transition. The Yang-Lee zeros at 9 = ±9 e are called the 
Yang-Lee edge zeros whose locations have never been known exactly for T > T c and d > 2. 
At the critical temperature the gap disappears, i.e., 9 e = 0. As temperature changes from 
T c to oo, the Yang-Lee edge 9 e moves from 9 = to 9 = n 

In addition to Yang-Lee zeros, Fisher showed that the thermal properties of the square- 
lattice Ising model for H = are completely determined by the zeros of the canonical 
partition function Z(T) in the complex temperature plane (Fisher zeros) [3, S|- Until now, 
Yang-Lee zeros and Fisher zeros (and other kinds of zeros for some problems) of numerous 
physical, chemical, biological, and mathematical systems have been investigated extensively 
to understand their important properties 



If the density of zeros 
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17| is found, the free energy, 



the equation of state, and all other thermodynamic functions can be obtained. However, very 
little is known about the actual form of the density of zeros. The density of Yang-Lee zeros 
: or the Ising ferromagnet has never been known exactly except for the one-dimensional case 
3]. Suzuki et al. Q| obtained the exact grand partition functions of the Ising model on the 
4x6 square lattice and the 3x3x3 simple-cubic lattice, and calculated the densities of Yang- 
Lee zeros for these small-size lattices. Kortman and Griffiths [3] investigated the densities of 
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Yang-Lee zeros for the Ising ferromagnet on the square lattice and the diamond lattice, based 

on the high-field, high-temperature series expansion, and found that the density of zeros for 

the square-lattice Ising ferromagnet diverges at Yang-Lee edge 9 e in high temperatures. For 

n 

the LxL square-lattice Ising model, Creswick and Kim [4J obtained the exact grand partition 
functions up to L = 10 and the semi-exact grand partition functions up to L = 14, calculated 
the densities of Yang-Lee zeros for T < T c , and evaluated the spontaneous magnetization 
TOo(T) = 1-ng(Q = 0,T) in the infinite-size limit. From Monte Carlo data, the density of 
Yang-Lee zeros near the transition point has been studied for the four-dimensional Ising 
fe mmag „e t Q a nd the lattice XY fe_t HQ 

. The density of Yang-Lee 



zeros has also been investigated experimentally for the two-dimensional Ising ferromagnet 
FeCl 2 in axial magnetic fields Q ■ 

Until now, the divergence of the density of zeros at Yang-Lee edge (the so-called Yang- 
Lee edge singularity) for the square-lattice Ising ferromagnet in high temperatures has been 
obtained only by the series expansion , 19]. Furthermore, for T > T c , the finite-size 

effects of the density of Yang-Lee zeros for the Ising ferromagnet have never been studied. In 
this paper, we investigate the properties of the density of Yang-Lee zeros based on the exact 
grand partition functions of the LxL square-lattice Ising model (L = 3 ~ 16). We deal 
with the three different classes of phase transitions for the Ising ferromagnet, (1) first-order 
phase transition, (2) second-order phase transition, and (3) Yang-Lee edge singularity, in 
a unified framework using the densities of zeros for finite-size systems. We also study the 
density of Yang-Lee zeros and its divergence in high temperatures from the finite-size data. 

II. DENSITY OF ZEROS 

For a finite-size system of size L, the density of Yang-Lee zeros g(6,L) at the critical 
temperature T c is given by P , 

g(6,L) = L~ d+yh g(6L yh ) (1) 

near 6 = 0. Equation (1) implies 

g (6) ~ \e\( d -y^ = {el 1 / 6 (2) 

for the infinite system. At the critical temperature the value of the magnetic scaling exponent 
Dh for the Ising ferromagnet in two dimensions (d = 2) is well known to be 15/8, that is, 



5 = 15. Therefore, the density of zeros vanishes at the critical point H c = (6 = 0). This 
kind of behavior of the density of zeros at T c is the characteristic of a second-order phase 
transition. On the other hand, for T < T c , we have the finite density of zeros at 9 = 0, 
g(0) = mo/2n, which is the characteristic of a first-order phase transition. 

To describe the density of zeros above the critical temperature as well, we propose that 
Eqs. (1) and (2) can be generalized as 

g{Q,L) = L- d ^g{Qm) (3) 

and 

g{Q) ~ \Q\( d -y^/y\ (4) 

where 6 = 9 — 6 e . For T < T c , Eq. (3) is the same as Eq. (1) because of 9 e = 0. The 
form g(Q) ~ |0| <T was suggested by Suzuki |l7[ who assumed o = 1/5 = for the two- 
dimensional Ising ferromagnet even above the critical temperature. According to Suzuki's 
assumption, the density of zeros above the critical temperature vanishes at Yang-Lee edge 

e e . 

Using the high-temperature series expansion for the Ising ferromagnet on the square lat- 
tice, Kortman and Griffiths [3] showed that the density of Yang-Lee zeros, g(Q) ~ 1 



diverges with a = —0.12 
tion. Kurtze and Fisher 



_5) at 6 = 6 e for T > 2T C , in disagreement with Suzuki's assump- 
18| refined the estimation of a by analyzing the high-temperature 
series expansion for the classical n-vector model and the quantum Heisenberg model in the 
limit of infinite temperature, and reported a = —0.163(3) in two dimensions (square and 
triangular lattices). Baker et al. analyzed the series expansion of much greater length 
for the square-lattice Ising ferromagnet, and obtained o = —0.1560(33) at y = e~ 2/3J = | 
(T w bJ/ks) and a = —0.1576(34) at y = g (T QJ/ks) from integral approximants. 
Remarkably, the exponent a has also been experimentally measured to be —0.15(2) in the 
range 49K < T < 53K and —0.365 at T = 34i^ for the triangular-lattice Ising ferromagnet 
FeCl 2 Q- 

Fisher renamed the edge zero as the Yang-Lee edge singularity for T > T c , and 
proposed the idea that the Yang-Lee edge singularity can be thought of as a new second- 
order phase transition with associated critical exponents. Above the critical temperature 
the value of a is known to be a = — | in two dimensions j^jj, resulting in yh = y (according 
to Eq. (4)), clearly different from the value of yn — y for the Ising ferromagnet at the 



critical temperature. The study of the Yang-Lee edge singularity has been extended to 
the spherical mode l |22l| . the quantum one-dimensional transverse Ising model 



hierarchical model [2J], branched polymers 
;he Ising systems with correlated disorder 
281 ] . and the antiferromagnetic Ising model 



the 

25l |. the Ising models on fractal lattices j^j], 




fluid models with repulsive-core interactions 
etc. 



III. GRAND PARTITION FUNCTION AND NUMBER OF STATES 

The Ising model in an external magnetic field if on a lattice with N s sites and iVj, bonds 
is defined by the Hamiltonian 

H = J E(l - Wi) +^E(i- (5) 

where J is the coupling constant, indicates a sum over all nearest-neighbor pairs of 

lattice sites, and o\ = ±1. The grand partition function of the Ising model is 

W„} 

where {cr n } denotes a sum over 2 Ns possible spin configurations and j3 = (/cbT) -1 • If we 
define the number of states, Q(E, M), with a given energy 

E = \Y,^-^3) (6) 

{hi) 

and a given magnetization 

M = V(1-^), (7) 
z i 

where E and M are positive integers < E < and < M < N s , then the grand partition 
function can be written as 

Z(x,y) = j: Y,n(E,M)y E x M , (8) 

E=Q M=0 

where y = e~ 2f3J and x = e~ 2f3H . From Eq. (8) it is clear that Z(x, y) is simply a polynomial 
in x and y. 

The exact integer values for the number of states Q(E, M) of the Ising model on the Lx L 
square lattices with cylindrical boundary conditions have been evaluated up to L = 10 to 
study the density of Yang-Lee zeros for T < T c using the microcanonical transfer matrix 



(/iTM) |4|. In Ref. [4], for 11 < L < 14, we had to evaluate the coefficients, for a fixed value 
of y, 

u(M) = Y / ^(E,M)y E 

E 

as real numbers of finite precision due to memory limitations. This semi-exact method 
worked well for T < T c , and produced all zeros on the unit circle. However, for T > T c 
the number of zeros on the unit circle, obtained from Z(x) = J2^(M)x M , decreases as 
T increases. Therefore, to study the density of Yang-Lee zeros for T > T c we need the 
exact integer values for Q(E,M), not u(M). In this work, the exact integer values for 
Q(E,M) are calculated up to L = 16 (2 256 = 1.158 x 10 77 configurations). The memory 
requirement during /xTM calculation is 2 x 4 x (N s /2 + 1) x (N b + 1) x P rn x 2 L byte for the 
two-dimensional Ising model, where P m is the maximum size for storing very long integer 
numbers. For L = 16, the memory requirement is 3.03 x 10 11 byte. (For comparison, the 
memory requirement for L = 10, the largest size calculated in Ref. [4], is 3.19 x 10 s byte. The 
largest size calculated before this work is L = 14 whose memory requirement is 3.44 x 10 10 
byte |30j.) The largest number of states for L = 16 is 

fi(249, 128) = 207964931539339552089271840638951892633 

537987437849943362668628644700590272, (9) 

which is approximately 2.08 x 10 74 . 



IV. NUMERICAL RESULTS 

For a finite-system of size L, the density of zeros (per site) may be defined as 

9 (M^^v (10) 

where 9k (k = 1,2,3,...) is the argument of the k-th zero on the unit circle Xq = e td of 
Z(x,y). The index k is counted from x = 1 (9 = 0) to x = — 1 (9 = n) for the zeros 
with 9 > 0. Hence, 9\{L) = 9 e (L). Figure 1 shows the density of Yang-Lee zeros for 
the 16 x 16 Ising ferromagnet at various temperatures (in unit of J/ks)- The Yang-Lee 
zeros at T = are uniformly distributed on the unit circle |6], and the density of zeros 
is given by g(9) = l/2ir (= 0.159...). As shown in Fig. 1, for a low temperature T — 1, 
the density of zeros is nearly identical with g(9) = 0.159. The spontaneous magnetization 



can be evaluated from the finite-size data for the density of zeros (mo(L) = 27rg(0,L), 



L = 3 ~ 16) using the BST extrapolation algorithm 3l(. For T = 1 the extrapolated value 
of the spontaneous magnetization is m = 0.9992757(3) in excellent agreement with the 
exact result 0.9992757519... j^j. At the critical temperature T c (= 2.269...), the maximum 
value of the density of zeros is 0.166096 at 9 = 2.541 where the density of zeros begins to 
decrease monotonically as 9 decreases. At T c the density of zeros vanishes at 9 = in the 
limit L — > 00, according to Eq. (2). For T = 3 the density of zeros is similar to that at T c , 
except for a violation of the monotonic decrease of the density of zeros around 9 = 0.2. It is 
interesting that the minimum value of the density of zeros at T = 3 (g(9 e = 0.126) = 0.112) 
is lower than that at T c ((7(0.023) = 0.117). Therefore, we may guess that the density of 
zeros for T just above T c , for example T = 3, vanishes at Yang-Lee edge 9 e in the infinite-size 
limit. As shown in Fig. 1, at high temperatures such as T = 5 and 10, the overall density 
of zeros becomes higher because the Yang-Lee edge 9 e tends to 7r as T rises. For T = 5 the 
density of zeros has a local maximum g = 0.182 near 9 e = 0.532. Finally, for T = 10, the 
density of zeros has the global maximum g = 0.279 near 9 e = 1.161. This fact may be a 
precursor for the divergence of the density of zeros at 9 e in the infinite-size limit for high 
temperatures. 

Figure 2 shows the densities of Yang-Lee zeros for L = 12 and 16 at different temperatures. 
Well below the critical temperature (for example, T = 1), the densities of zeros for L — 12 
and 16 are not distinguishable, as shown in Fig. 2(a). Therefore, the finite-size effect for 
the density of zeros is very small well below T c . At and near the critical temperature 
(for example, T = 2.5 in Fig. 2(a)), the densities of zeros for L — 12 and 16 are also 
indistinguishable, except for them near 9 = where the density of zeros g{L) decreases 
sharply as the system size L increases. Well above the critical temperature (for example, 
T = 10 in Fig. 2(b)), the overall forms of the density of zeros for L = 12 and 16 are 
almost identical. The finite-size effect for the density of zeros at T = 10 is relatively large, 
compared to those at T = 1 and 2.5. Near the Yang-Lee edge 9 e = 1.161, the density of 
zeros is somewhat irregular, and it increases as L increases. 

If we set 6 = (0 = 9 e (L)) in Eq. (3), we have 

g(0 1 L) = L^g(0) } (11) 
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from which we can define the magnetic scaling exponent 

ln[g(0,L + l)/g(0,L)] 
^ (L) " d+ kp + 1)/L] (12) 

for finite lattices. Table I shows the values of yh{L) at T = 1, T c , and 10, respectively. For 
T = 1, Dh{L) changes very slightly around d = 2 as L increases. This kind of behavior 
may result from the small finite-size effect for the density of zeros well below the critical 
temperature. At T c and T = 10, yh{L) increases monotonically as L increases. In particular, 
at T = 10, yh(L) increases quickly. 

Table II shows the extrapolated results for the magnetic scaling exponent by the BST 
algorithm. For T — 1, the extrapolated result clearly indicates yh — d — 2 that is expected 
at a strong first-order phase transition j^. At T c , the extrapolated value of yh = 1.8747(7) 
is in agreement with the exact value of yh = y — 1.875. For T > 3, we obtained the 
extrapolated value of yh = 2.4 which gives a = — ^ (the characteristic of Yang-Lee edge 
singularity). According to Eq. (4), the density of zeros diverges at Yang-Lee edge 9 e when 
yh > d. Finally, we can easily distinguish the three different classes of phase transitions 
(first-order phase transition, second-order phase transition, and Yang-Lee edge singularity) 
by estimating yh from the densities of zeros on finite lattices. 

Near the critical temperature T c = 2.269..., we expect some crossovers between different 
classes. At T = 2.2 we obtained yh = 1.98(15) that is close to yh — 2 but is not clearly 
distinguished from yh = -j due to the large error. At T = 2.4 and 2.5, the extrapolated 
values for y h are 2.33(23) and 2.89(72), respectively, indicating the existence of Yang-Lee 
edge singularity. But their errors are too large. It is interesting that the densities of zeros 
for finite-size systems at T = 2.4, 2.5, and 3 are quite similar to those at T c but they give 
the values of yh (> d = 2) completely different from yh = y (< d). For the first time, it 
is shown that, even in the range of 2.4 < T < 4, the square-lattice Ising ferromagnet has 
the Yang-Lee edge singularity. At T = 2.3 we obtained yh = 1.58(4) whose nearest value is 
Vh — y a ^ This value of yh = 1.58(4) implies that the density of zeros for T > T c vanishes 
at 9 e . However, if the data of much larger sizes can be used, there will be a possibility of 
obtaining yh > 2 even at T = 2.3. 

The finite-size data for the density of zeros g(9 e ,L) (L = 3 ~ 16) in high temperatures 
have also been extrapolated to infinite size. In high temperatures, it is convenient to deal 
with the inverse density of zeros f(9 e , L) = l/g(9 e , L) because f(9 e , L) vanishes with the scale 
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factor L~( yh ~ d \ in the limit L — > oo, but g(9 e , L) diverges. Table III shows the extrapolated 
results for the inverse density of zeros f(9 e ). For T > 5 the extrapolated results imply the 
divergence of the density of zeros at 9 e . However, at T = 4, we can not conclude anything 
only from the extrapolated result for the inverse density of zeros because of the large error. 



V. PHASE TRANSITIONS IN SMALL SYSTEMS 



In the previous section, the density of zeros has been extracted from the locations of Yang- 
Lee zeros for finite sizes, using Eq. (10) P. liol [lfl . Recently, two alternative approaches (JK 



approach 



13 



35 



361 ] and BMH approach |36j, |37() have been proposed to extract the 



density of zeros from the locations of Fisher zeros in the complex temperature (y = e~ 2 @ ) 
plane for finite sizes and to study phase transitions in small systems. These two approaches 
are based on the assumption that Fisher zeros i/j = Sj + itj = y c + irj exp(-iip) close to the 
critical point y c lie on a single line for quite large L, where the index j(= 1, 2, 3, ...) increases 
with distance rj from the critical point, Sj = Ke(yj) = y c + rjsmip, and tj = lm(yj) = 

Tj COSIp. 

In the JK approach 0,0, HQ, the density of zeros is defined as 

X L (r) = L- d j:S(r-r 3 (L)) J (13) 

i 

and the cumulative density of zeros A^(r) = J r \L{r')dr' becomes a step function with 



Arlr) 



J_ 

L d 



[if r e {rj,r j+1 )} 



(14) 



for a finite-size system of size L. Then the average cumulative density of zeros is given by 

2J-1 



2L d 



and it is expressed as 



Aoo(^) = Aoo(0)r + ar 



w+l 



+ 



(15) 



(16) 



in the limit L — > oo. Here, the density of zeros A oo (0) on the real axis is proportional to the 
latent heat Ae. At a first-order phase transition, the first term Aoo(0)r becomes dominant 
for small r. At a second-order phase transition, the second term ar w+l is the leading one 
because of Aoo(0) = 0. Finally, the cumulative density of zeros can be approximated as 



2L d 



Al(tv) = ai[t j (L)] a2 + a 3 



(17) 



because the parameter rj for Fisher zeros close to the critical point y c may be taken to be the 
imaginary part tj(L) due to y c ~ Sj(L). A value of inconsistent with zero indicates the 
absence of a phase transition. The parameter a 2 determines the order of a phase transition. 
A first-order phase transition takes a value a 2 ~ 1 for small t, whereas a value of a 2 larger 
than 1 indicates a second-order phase transition whose specific-heat exponent is given by 
a = 2 — a 2 . For a first-order phase transition, the parameter ai corresponds to the density of 
zeros Aoo(0) on the real axis. The parameter a 2 is also expressed as a 2 = d/yt in terms of the 
thermal scaling exponent y t . The JK app roach has also been applied to the Yang-Lee zeros 
of the square-lattice XY ferromganet where the parameter a 2 is expressed as a 2 = d/yh 
in terms of the magnetic scaling exponent y^, that is, 

2d , . 

a 2 = . 18 

d+2-r] y ' 

The JK approach has also been generalized to study the areal distributions of Fisher zeros 

On the other hand, there are three important parameters in the BMH approach 0, S| • 
The first parameter t\ (the imaginary part of the first zero) becomes in the limit L — > oo 



for a true phase transition in the Ehrenfest sense 3, y, q, y, Q- 

If there is a ph 

transition, we assume that the line density of zeros, 



asc 



Mi) = ~ ( I " 1 + 1 " r ) , (19) 

2 VI%-%-il \yj+i-yj\J 

follows a simple power law 4>(t) ~ t ^ for small t. The second parameter can be estimated 
as 

In0(t 3 )-ln0(t 2 ) 

a * = lnt 3 -lnt 2 (20) 

from the first four Fisher zeros (yi, V2, Vs, Vi) ■ The line of Fisher zeros also crosses the real 
axis at the critical point y c with the angle ip, yielding the third parameter 

7 V , = tan?/> = ^ — Ji. (21) 
£2 — £1 

The second and third parameters (a^ and 7^,) determine the order of a phase transition. The 
values of a<f> = (a^ can be less than zero for finite sizes) and 7^ = indicate a first-order 
phase transition. A second-order phase transition occurs for < a$ < 1 and arbitrary 7^, 
and a higher-order phase transition for > 1 and arbitrary 7^. Alves et al. 3(| have 
performed numerical comparisons of the JK and BMH approaches for the Fisher zeros of 
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the two-dimensional four-state and five-state Potts models and of biological molecules, using 
Monte Carlo data. Their results implied that the JK approach is little better than the BMH 
approach. 

The approach explained in the previous sections, which may be called as the angular- 
density-of-zeros (ADOZ) approach, can also be used to study phase transitions in small 
systems. In Eq. (4), the density of Yang-Lee zeros near = is given by g(Q) ~ \Q\( d -yh)/vh _ 
Similarly, the density of Fisher zeros near = is also given by g(Q) ~ \Q\( d -y*)/yt [16]. 
The exponent (d — yt)/yt (or (d — yh)/yh) corresponds to the parameter a 2 — 1 of the JK 
approach. As in the BMH approach, the zero value of t\ in the limit L — > oo indicates the 
existence of a true phase transition. When there is no phase transition, Fisher (or Yang- 



Lee) edge singularity [39|, |4Cj exists if y t > d (or y h > d). As shown in Table I, the values 



of yh{L) at T = 10 are clearly greater than d = 2 for relatively large L, suggesting the 
existence of Yang-Lee edge singularity. When there is a phase transition, the value of y t 
(or yh) determines the order of the phase transition. To calculate yt(L) (or yh(L)), we need 
only the first two Fisher (or Yang-Lee) zeros according to Eqs. (10) and (12). A value of 
yt ~ d (or yh ~ d) means a first-order phase transition. The values of yh{L) at T = 1 
in Table I clearly indicate a first-order phase transition. A second-order phase transition 
occurs for d/2 < y t < d (or d/2 < yh < d), and a higher-order phase transition for y t < d/2 
(or y h < d/2). The values of yu{L) at T c in Table I indicate a second-order phase transition. 

Now the three (JK, BMH, and ADOZ) approaches are applied to Yang-Lee zeros of the 
square-lattice Ising ferromagnet. The BMH approach is applied to Yang-Lee zeros for the 
first time. Because the JK and BMH approaches assume a line distribution of zeros near 
a transition point, we use Yang-Lee zeros from L = 8. The JK and BMH approaches have 
been compared for the approximate zeros of quite large (but still small) finite-size systems 
using Monte Carlo simulation [3iJ. Here these approaches are compared for the exact zeros 
of very small finite-size systems {L — 8 to 16). Because the BMH approach requires the first 
four zeros, we use only them in the applications of the JK and BMH approaches for a fair 
comparison. The ADOZ approach requires the first two zeros of two different sizes, in total, 
four zeros (see Eq. (12)). 

In the identification of the orders of phase transitions in small systems, the most im- 
portant quantities are a 2 (L), a^(L), and yh(L), respectively, for the JK, BMH, and ADOZ 
approaches. Figure 3 shows the results of cli{L), a^(L), and yh{L) at T = 1 (first-order 
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phase transition), T = 2.2 (weak first-order phase transition), and T c = 2.269... (second- 
order phase transition) for the Yang-Lee zeros of the square-lattice Ising ferromagnet. Here, 
a weak first-order phase transition means a first-order phase transition with very large (but 
finite) correlation length. All the approaches clearly identify the first-order and second- 
order phase transitions. In particular, the finite-size results at T = 1 are perfect in that 
they already reach the infinite-size values a 2 = 1, = 0, and yh = 2. But the three 
approaches fail to identify the weak first-order phase transition. An advantage of the JK 
approach is that we can use the zeros from different sizes together in the estimation of a 2 . 
The only first Yang-Lee zeros from L — 13 to 16 (in total, four zeros) are used to estimate 
a 2 at T = 2.2 and T c . At T c , we obtain a 2 = 1.0874(4), close to the infinite-size value 
a 2 = d/y h = 16/15(= 1.0666...). At T = 2.2, the result is a 2 = 1.011(2). It is remarkable 
because the JK approach almost identifies the weak first-order phase transition only from 
very small systems up to L = 16. 



VI. CONCLUSION 



We have enumerated the exact integer values for the number of states Q(E, M) of the 
Ising model on the L x L square lattice up to L = 16, from which the grand partition 
function Z(T,H) = Ei?,M fi (- E . M )e" 2/3(J£+ ™ ) is obtained and the density of Yang-Lee 
zeros g(6, L) on the unit circle x = e~ 2l3H = e l6 is evaluated for the Ising ferromagnet. We 
have compared the densities of Yang-Lee zeros at different temperatures for L = 16. We 
have also discussed the finite-size effects for the densities of zeros at different temperatures. 
The three different classes of phase transitions for the Ising ferromagnet, (1) first-order 
phase transition, (2) second-order phase transition, and (3) Yang-Lee edge singularity (the 
divergence of the density of zeros at Yang-Lee edge in high temperatures) have been clearly 
distinguished by estimating the magnetic scaling exponent yh from the densities of zeros 
on finite lattices. The Yang-Lee edge singularity that has been detected only by the series 
expansion until now for the square- lattice Ising ferromagnet, has also been obtained using 
the densities of zeros for finite-size systems. We have also discussed the identification of 
the orders of phase transitions in small systems using the three (JK, BMH, and ADOZ) 
approaches to extract the density of zeros from numerical data. 
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TABLE I: The magnetic scaling exponent yh{L) for finite lattices, obtained from the density of 
zeros g(6 e ,L) using Eq. (12), at T = 1, T c , and 10. 



L 


T = 1 


T c 


T = 10 


3 


1.994240 


1.776068 


1.676531 


4 


1.999877 


1.795022 


1.819092 


5 


2.001264 


1.806392 


1.917883 


6 


2.001516 


1.814538 


1.992971 


7 


2.001475 


1.820774 


2.051768 


8 


2.001368 


1.825750 


2.098558 


9 


2.001256 


1.829834 


2.136280 


10 


2.001153 


1.833257 


2.167057 


11 


2.001062 


1.836174 


2.192455 


12 


2.000982 


1.838693 


2.213638 


13 


2.000913 


1.840891 


2.231483 


14 


2.000852 


1.842827 


2.246655 


15 


2.000798 


1.844547 


2.259666 
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TABLE II: The extrapolated results for the magnetic scaling exponent at different temperatures. 



T 


iii, 

yn 


1 


2 000000(6) 


2.2 


1 98(15) 


T 


1.8747(7) 


2.3 


1.58(4) 


2.4 


2 33(23) 


2.5 


2 89(72) 




2 3Q(3) 


4 


2.4056(11) 


5 


2.4005(3) 


10 


2.40026(5) 


15 


2.4001(2) 


20 


2.3998(3) 



TABLE III: The extrapolated results for the inverse density of zeros f(0 e ) at Yang-Lee edge 9 e in 
high temperatures. 



T 6? f(6 e ) 

4 0.314846(75) 0.27(756) 

5 0.514241(33) -0.01(7) 
10 1.149334(14) -0.05(4) 
15 1.479519(10) -0.02(4) 
20 1.687041(8) -0.005(31) 
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FIG. 1: (Color online) The density of Yang-Lee zeros for the 16 x 16 Ising ferromagnet at T = 1, 
T c (= 2.269...), 3, 5, and 10 (in unit of J/k B ). 
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FIG. 2: (Color online) The density of Yang-Lee zeros for L = 12 and 16. (a) T = 1 and T = 2.5. 
(b) T = 10 (in unit of J/fc B )- 

18 



1.25 



1.20 



1.15 



a 2 (L) 
1.10 



1.05 - 



1.00 




0.15 


- (b) 










0.10 






* 










T=Tc 
T=2.2 








0.05 


■ 


T=1 








0.00 














8 


10 


12 
L 


14 


16 


2.05 












: (c) 










2.00 












1.95 


- • 


T=Tc 








y h (L) 


■ 


C\J 

II II 








1.90 












1.85 






• * 


• • 


• 



10 



12 



14 



16 



FIG. 3: (Color online) The values of (a) ai (L) (JK approach), (b) a^L) (BMH approach), and 
(c) yh(L) (ADOZ approach), respectively, for T = 1, 2.2, and T c (in unit of J/ks), obtained from 
Yang-Lee zeros of sizes L = 8 to 16. 
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